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ABSTRACT 


Thirty-six  XBT  temperature  profiles  have  been  used  in  a  parametric  model  introduced  by  Hendry 
(1988)  to  model  die  Gulf  Stream's  thermal  structure  at  65"W  between  200  and  1200  dbar,  with  an  ms 
residual  error  of  0.56°C.  Velocity  has  been  computed  geostropbically  relative  to  1200  dbar,  and  has 
been  included  in  calculating  potential  vorticity  analytically  from  the  model.  The  resulting  potential  vor¬ 
ticity  section  for  65"  W  has  been  compared  with  the  analogous  result  from  Hendry’s  parametric  model 
at  59" W,  as  well  as  observed  potential  vorticity  sections  from  68°  to  55°  W.  There  is  a  significant 
feature  in  the  potential  vorticity  structure  at  65°  W  not  found  at  59°  W,  namely  a  relative  minimum  in 
potential  vorticity  along  isopycnals,  catered  at  the  Gulf  Stream's  axis  and  350  dbar.  The  modelled 
potential  vorticity  sections  are  consistent  wife  fee  observations,  including  fee  downstream  disappearance 
of  this  feature.  The  dynamical  implications  of  these  results  are  briefly  discussed. 
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1.  Introduction 


Various  authors  (Hall  and  Bryden,  1983;  Hendry,  1988;  Hogg,  1986,  1991)  have  used  parameteri- 
rjirinn*  of  die  Gulf  Stream’s  (GS)  thermocline  temperature  field,  based  on  limited  measurements,  in 
order  to  correct  current  meter  velocity  and  temperature  time  series  to  fixed  standard  depths.  These 
parameterization  generally  take  advantage  of  the  quasi-fixed  or  "canonical"  nature  of  the  field,  at  least 
for  a  given  downstream  location:  the  primary  mode  of  variability  for  the  GS  thermocline  is  a  simple 
horizontal  translation,  rather  than  internal  deformation  (Manning  and  Watts,  1987;  Halkin  and  Rossby, 
1983;  Rossby,  1987).  Beyond  their  utility  as  correctional  tools  these  parameterizations  produce  model 
GS  thennoclines  that  can  themselves  be  analyzed.  Indeed,  in  his  1988  paper,  Hendry  suggests  using 
them  for  regional  intercomparisons  of  GS  structure.  The  intent  of  this  paper  is  to  do  precisely  that 
using  data  collected  in  1988. 

During  a  hydrographic  survey  made  as  part  of  the  Synoptic  Ocean  Prediction  (SYNOP)  Experi¬ 
ment  (Hall  and  Fofonoff,  1993)  four  XBT  crossings  of  the  GS  were  obtained,  between  additional  CTD 
and  XBT  sections  at  68°  and  33°W.  Hall  and  Fofonoff  (1993,  hereafter  HF93)  have  described  in  detail 
the  observed  GS  structure  from  foe  CTD  sections,  but  they  did  not  discuss  die  XBT  sections.  Three  of 
these,  located  from  67°-63°W,  consisted  entirely  of  T5  drops,  which  yield  temperatures  as  deep  as 
1800  m;  they  were  acquired  within  a  time  span  of  just  over  a  day,  and  crossed  die  "straight"  portions  of 
a  trough/crest  pair  (Fig.  la).  This  group  of  sections,  comprising  a  total  of  36  temperature  profiles,  has 
adequate  vertical  and  cross-stream  resolution  to  apply  Hendry’s  model  and  obtain  a  smoothed  GS  ther¬ 
mocline  nominally  at  65°  W,  for  comparison  with  his  39°W  section. 

In  order  to  make  this  comparison  with  Hendry’s  GS  thermocline  at  39°W,  Hendry’s  (1988) 
method  is  applied  to  foe  set  of  XBT  profiles  centered  around  63°W.  La  the  next  section,  the  model  fit 
is  derived  in  some  detail,  and  the  resulting  analytical  GS  thermocline  is  compared  with  foe  XBT  obser¬ 
vations  and  with  Hendry's  59°  W  section.  In  section  3,  potential  vortidty  cross-sections  for  this  65°  W 
and  Hendry’s  39®  W  section  are  computed  to  intercomparison.  In  this  paper,  potential  vortidty  will  be 
foe  dynamical  quantity  of  particular  interest,  to  it  is  what  determines  the  current’s  inherent  stability 
properties.  Does  the  OS's  canonical  cross-stream  structure  change  in  potentially  significant  ways  as  we 
move  downstream?  The  discussion  places  foe  results  into  a  larger  context  by  comparison  with  observa¬ 
tions  of  GS  potential  vortidty  at  other  locations  along  foe  Stream  as  well.  The  final  sections  briefly 
summarize  the  conclusions  of  foe  overall  study. 


2.  TfcaMnM 

A  brief  overview  of  Hendry's  method,  roughly  following  his  notarial,  is  is  follows:  Hist,  each 
temperature  profile  was  fit  to  a  hyperbolic  tangent  function  of  pressure  over  a  range  of  200-1300  dbar 

T  -  A  tanh[B(C-p)]  +  D  (1) 

so  that  (A  tanh(BCy*D)  and  (D-A)  bound  the  temperatures  a  given  profile  can  represent  (since  p  >  0 
only);  C  is  die  "thermociine  depth,"  where  dT/dp  is  a  maximum;  and  B  gives  the  scale  depth  of  the 
tbermocline.  Next,  parameters  A,  B,  and  D  were  fit  to  functions  of  C,  viz. : 


i  •  ye(rtC)  +  e; 

(2) 

D  m  aA  +  P; 

(3) 

a  sech[b(C-c )]  +  d. 

(4) 

and  finally,  with  cron-stream  distance  y  taken  as  the  distance  to  the  location  of  the  15° C  isotherm  at 
200  m,  parameter  C  was  fit  to 

C(y )  -  -Aj  tanhVtiCry)]  +  D,  (5) 

where  the  subscript  T’s  are  used  to  indecate  that  these  are  different  from  (A,  B,  C,  D)  in  (1).  In  (3), 
(Ai  +  Dj)  and  (Dl  -  At)  give  the  maximum  and  minimum  themocline  depths;  flj  is  foe  scale  width  of 
the  thermociine  drop,  U. ,  the  current’s  deformation  radius;  and  C,  is  an  offset  dependent  on  origin 
choice.  Positive  y  values  correspond  to  increasing  distance  onshore  (northward)  of  foe  axis.  Thus,  x 
and  u,  foe  alongstream  coordinate  and  velocity,  are  positive  in  in  foe  downstream  direction,  usually 
date  to  east;  and  v  and  y,  the  cross-stream  coordinate  and  velocity,  are  positive  in  die  "northward-like” 
directions,  onshore.  (In  most  of  foe  literature  on  foe  Gulf  Stream,  foe  onshore  side  is  shown  to  foe  left 
in  figures  of  cross-seed ons,  so  this  convention  is  adopted  for  the  figures  in  this  paper.)  The  notation  is 
admittedly  confusing,  but  reference  to  Table  1  shows  which  of  these  parameters  are  constant  rather  than 
fimerions  of  y.  Previous  authors  (for  example,  Hogg,  198 6)  have  used  foe  hyperbolic  tangent  T(p) 
profile  locally  as  a  GS  model;  but  the  temperature  range  and  scale  depth  of  foe  GS  thermociine  change 
•croaa  the  current,  and  fins  cron  stream  structure  must  be  included  to  obtain  a  realistic  representation. 
[Clearly,  it  might  be  desirable  to-  find  a  model  depending  on  fewer  than  13  parameters;  however,  the 
point  of  fins  note  is  instead  to  provide  another  cross-section  to  compare  with  Hendry’s.] 

Figures  l(byc)  show  the  collection  of  temperature  profiles  from  the  36  XBT  drops  to  be  used  in 
die  model,  rod  fin  temperature  section  verms  along-track  distance  from  southwest  (right,  0  km)  to 
northeast  (left,  300  km).  There  are  several  temperature  inversions  associated  with  Slope  Water 
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intrusions  on  the  cold  "sides"  as  well  as  a  deep,  thickened  18 °C  Water  layer.  These  features  will  con¬ 
tribute  to  the  uncertainty  in  the  final  fits,  but  there  is  no  good  way  to  quantify  their  specific  effect 
Errors  in  fail  rate  for  TS’s  are  another  source  of  uncertainty,  because  these  errors  have  a  greater  relative 
effect  deeper  in  die  water  column,  where  the  net  change  of  temperature  is  small  across  the  GS  at  a 
given  depth,  the  following  fit  includes  values  only  between  200  and  1200  dbar,  rather  than  1300  dbar  as 
Hendry  used.  Proceeding  directly  along  Hendry’s  lines,  each  of  the  36  temperature  profiles  was  then  fit 
to  equation  (1),  yielding  a  collection  of  36  values  of  (A,  B,  C,  D)  corresponding  to  the  36  XBT  profiles. 
(Meters  nod  decibars  assumed  interchangeable  for  the  purposes  of  this  note.)  An  initial  fit  of  A  to  C 
was  used  to  identify  values  lying  more  than  1.3  standard  deviations  from  the  fitted  curve.  This  elim¬ 
inated  five  points,  and  die  remaining  31  stations  were  used  to  continue  with  the  fits  for  B(C)  and  D(A); 
these  fits  had  no  outliers  beyond  two  standard  deviations  (2.1  for  D).  Figures  2(a,  b,  d)  show  the  fitted 
curves  A(Q,  B(Q,  D(A);  the  31  points  to  which  the  curves  were  fit;  and  the  3  outlier  stations,  plotted 
as  open  circles. 

hi  order  to  calculate  cross-stream  position  y  as  the  distance  to  15°C  at  200  m,  one  must  define  y 
for  stations  lying  between  two  adjacent  stream  crossings  -  though  C  is  fixed  for  such  stations,  the 
choice  for  y  is  dearly  ambiguous.  [Note  that  using  a  different  definition,  such  as  12°C  at  500  m,  as  an 
origin  simply  shifts  die  y  coordinate  uniformly  by  about  17  km,  and  does  not  remove  the  ambiguities 
noted  above.]  Between  the  first  two  sections,  there  was  an  overlap  of  two  XBT  drops;  both  of  these 
yield  similar  (C,  y)  pairs,  and  both  were  included  in  the  fit.  However,  between  die  second  and  third 
sections,  we  concluded  on  the  basis  of  satellite  DR.  data  that  we  had  not  fully  crossed  die  GS  into  the 
warm  side  before  the  isotherms  rose  again  for  the  third  section.  There  are  four  stations  lying  in  the 
overlap  region  that  might  be  used  to  define  (C,  y)  pain  for  the  fit;  two  of  these  were  included  as 
"belonging  to"  die  second  section,  and  two  to  the  third.  These  choices  yield  a  group  of  33  (C,  y) 
values  for  the  fit  given  by  (3),  and  again  this  fit  was  good  to  within  two  standard  deviations.  Fig.  2c 
shows  die  fit  Ofy),  die  33  points  used  for  the  fit,  and  the  3  points  omitted  from  all  the  fits. 

The  mu  error  of  the  final  T(pj)  fit  to  die  data  for  all  36  stations,  between  200  and  1200  dbar,  is 
0-36*  C,  and  for  the  31  stations  used  in  the  fits,  it  is  0.30°  C.  The  values  of  the  thirteen  parameters  are 
listed  in  Table  1,  along  with  Hendry’s  (1988)  values,  and  Fig.  3  compares  die  synthesized  temperature 
field  with  die  first  XBT  temperature  section  over  a  similar  range  of  y.  Except  for  north  of  the  axis 
(y  >  0),  where  the  frontal  aspect  at  temperatures  T  £  14 °C  is  smoothed  by  the  model,  it  reproduces  the 
temperature  structure  well:  The  synthetic  field  has  a  tbennodine  depth  ranging  from  112  to  763  dbar 
and  reproduces  temperatures  ranging  from  4.27  to  18.98°  C  between  0  and  1200  dbar.  For  comparison, 
Hendry’s  model  tbermocline  varies  from  123  to  767  dbar  and  has  temperatures  ranging  from  4.32  to 


14.79 *C  for  0*1300  dbar.  The  similarity  is  hardly  surprising:  Recall  that  A  sod  D  define  the  net  tem¬ 
perature  range  of  die  current,  while  constants  Ax  and  D\  in  the  fit  for  C  give  the  net  tfaermodine  drop 
across  the  current,  and  B  t~‘  is  the  deformation  radius  for  the  thermocline,  approximately  38-40  km 
aocording  to  the  values  in  Table  1.  The  constant  parameter  values  in  these  three  fits  (Table  1)  are  only 
moderately  different  for  the  two  sections,  as  one  might  expect,  since  such  integral  properties  of  the 
current  are  relatively  constant  going  downstream. 

In  contrast,  the  local  inverse  scale  depth  for  temperature  is  modeled  by  fit  (4)  for  B,  which  is 
much  different  for  the  two  sections;  because  C  is  a  function  of  y,  horizontal  scales  of  variability  are 
implicit  in  B(C)  as  well,  even  though  C(y)  is  similar  for  the  two  sections.  Figure  4  shows  [B(C(y)rl 
explicitly  as  a  function  of  y  across  the  current  for  63*  and  59°  W:  not  only  are  the  extrema  (given  by 
d,  a+d  in  (4))  farther  apart  at  63°  than  at  39°  W,  but  the  cross-stream  distance  over  which  most  of  the 
change  occurs  is  slightly  less  at  tire  former  as  well,  leading  to  a  steeper  profile  for  0~l  overall.  As  will 
be  seen  in  Section  3,  the  stronger  character  of  the  increase  of  B  (decrease  of  B~l)  at  65°  W  leads  to  a 
significant  PV  feature  that  is  absent  at  39°W;  and  we  will  see  that  the  feature  is  robust  to  errors  in  the 
fit 

Geostrophic  velocities  may  be  calculated  using 

du  -go.  3T 

W--STW  (6) 

Oo  •  m"5  °C~l 

where  T  is  used  as  a  proxy  for  a*,  and  a,  is  the  proportionality  constant  Keeping  in  mind  the  resolu¬ 
tion  of  the  front  and  that  T  is  not  a  particularly  good  proxy  for  density  when  T  2  12°C  (especially  on 
the  cold  side),  the  resulting  velocity  field  nevertheless  reproduces  features  typical  of  GS  velocity  sec¬ 
tions,  such  as  exhibiting  a  southward  shift  with  depth  of  the  velocity  maximum.  Table  2  makes  a  quan¬ 
titative  comparison  of  maximum  velocities  unmxt  and  transports  in  depth  ranges,  for  the  model  and  the 
observations:  model  values  for  transport  are  intermediate  among  the  observations,  while  is  on  the 
low  side,  probably  because  it  occurs  within  the  front  Transports  for  Hendry's  model  at  39°W  are 
again  similar  to  results  at  65°  W,  because  these  depend  on  the  net  drop  of  isotherms  across  tire  current 
rather  than  local  variability. 

We  now  turn  to  the  PV  structure  of  die  current,  which  is  of  greater  dynamical  interest 
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3.  Potential  Vortidty 

An  analytic  model  is  useful  for  examining  a  highly  derived  quantity  such  as  the  potential  vomdty 
(PV)  structure  associated  with  the  temperature  field,  since  the  calculation  can  be  done  analytically  for 
die  most  pvt  As  noted  above,  temperature  has  been  used  as  a  proxy  for  density  o«  using  a  fixed  value 
oto  of  the  proportionality  constant;  thus,  potential  vorticity  n  can  be  calculated  using  a*  as  the  conser¬ 
vative  tracer  in  fire  expression  for  n,  for  comparison  with  HF93.  To  review,  the  full  PV  n  of  a  water 
parcel  is  then  given  as 

n  -  iL  (2D  +  ®^)  Vo9  (7) 

Po 


(e.g.,  Pedlosky,  1979),  where  2D  is  planetary  vortidty  and  is  the  (three-dimensional)  relative  vorh- 
dty  of  die  flow.  The  minus  sign  ensures  positive  values  of  n  in  the  northern  hemisphere  for  large 
scale  flow,  hi  quasi-geoatrophic  theory,  (7)  reduces  to 


n 


dv  du.  d°e 

T&'wsr 


(8) 


and  in  the  eastward  GS,  vx  «  Uy,  so  that  term  may  be  dropped.  However,  because  zsopycnals  depart 
significantly  from  being  horizontal,  a  non-quasigeostrophk  term  must  also  be  retained  from  the  dot  pro¬ 
duct  in  (7),  namely,  ft*  (scale  analysis  can  be  used  to  show  that  this  term  may  be  of  die  order 
Po  oz  ay 

of  ,  see  HP93).  The  expression  for  n  then  becomes,  as  in  HF93, 

Po  oy  oz 


it  _  -1  r*  _  _  1  ^ 

-  -  -3^}-3F  ~  " 


(9a) 


“Iff  -  to)3**  _  S3  s3**) 

P r  wsr 


or,  using  (6)  and  writing  (9a)  in  terms  of  temperature. 


_  a,  du  9T  ^<*odu9T 
"  pT^  "  -pTST^y 


(9b) 


..  du  dr  j<*o2  ar* 
PTV  “  ly}W  V 


One  may  add  a  barotropic  velocity  component  at  1200  dbar  if  desired,  but  leaving  it  out  in  subsequent 
calculations  does  not  affect  the  qualitative  structure  of  the  PV  fields.  Therefore  only  baroclinic  veloci¬ 
ties  relative  to  1200  dbar  are  included  in  what  follows,  hi  the  definition  (9)  of  n,  the  only  term  that  is 


not  easily  calculable  analytically  is  du/3y,  since  obtaining  the  gec*trophic  velocity  in  the  first  place 
involves  integration  with  depth;  it  is  far  simpler  (cod  adequately  accurate)  to  calculate  tins  term  by 
differencing. 


The  PV  fields  have  been  calculated  for  this  section  and  for  Hendry's,  using  the  parameter  values 
of  Table  1  and  equations  (6)  and  (9);  they  are  shown  in  Figures  5a  and  6a.  It  is  important  to  recognize 
that  a  density-conserving  water  parcel  in  the  current  locally  feels  the  PV  distribution  along  isopycnals. 
Therefore  several  representative  isotherms  have  been  superposed  on  the  PV  in  Figures  5  and  6,  to  give 


the  reader  an  idea  of  how  the  cross-stream  PV  structure,  as  viewed  by  water  parcels,  varies  from  one 
part  of  the  thermocline  to  another.  Several  other  comments  are  in  order  regarding  the  three  different 

components  contributing  to  n  in  (9).  First,  note  that  /,  u,,  and  j—  all  vary  across  the  Stream; 
d<fa  dog 

changes  in  the  product  /_j_  are  dominated  by  changes  in  _j_,  and  0  is  insignificant  compared  to  the 


relative  vortidty.  Thus  foe  term  proportional  to  / 


dae 
*3 T 


(or  fTj)  is  referred  to  as  foe  stretching  com¬ 


ponent  since  it  mainly  reflects  compression  and  stretching  of  isopycnals  (isotherms).  In  contrast, 
do« 

changes  in  are  dominated  by  the  change  in  fire  relative  vortidty  -u,  itself.  This  term  is 

referred  to  as  the  relative  vortidty  component,  and  it  reinforces  stretching  north  of  foe  axis,  the 
cyclonic  side,  but  reduces  it  an  foe  anticyclonic  side  where  ~Uy  is  negative.  The  third  term  will  be 
called  foe  non-quasigeostropbic  (non-QG)  term,  and  is  everywhere  negative.  It  is  strongest  just  north  of 
foe  velocity  maximum,  and  like  foe  velodty  field,  its  maximum  slopes  southward  with  depth.  The 
stretching  component  dominates  PV  structure  at  both  sections,  as  demonstrated  by  comparison  of  Fig¬ 
ures  Sa  and  6a  (total  PV)  with  5b  and  6b  (stretching  PV).  The  relative  and  non-QG  contributions  are  of 
comparable  magnitude  to  one  another,  as  much  as  25%  as  great  as  the  stretching  terms;  because  of  their 
characteristic  structures,  they  tend  to  caned  in  foe  cydcnic  portions  of  the  current  and  reinforce  one 


another  on  foe  anticydonic  side. 


The  GS  has  long  been  recognized,  mid  modeled,  as  a  strong  PV  front,  separating  less  strongly 
stratified,  low  PV  Sargasso  Sea  waters  from  more  stratified  higher  PV  slope  waters.  This  structure  is 
visible  but  not  prominent  in  Figures  5  and  6  because  foe  strangest  part  of  foe  front  lies  above  200  m, 
and  is  largely  due  to  foe  transition  to  slope  water  T/S  characteristics,  not  accounted  for  by  foe  model 
A  more  complex  structure  is  striking  at  65"  W  (Fig.  5),  namely  a  pronounced  minimum  along  isotherms 
loomed  near  foe  axis  at  a  depth  of  about  350  m.  Notice  that  this  minimum  is  inherent  in  foe  vertical 
density  structure  alone  (Fig.  5b),  for  its  appearance  does  not  depend  on  foe  inclusion  of  foe  relative  and 
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nco-quasigeostrophic  terms  of  IT,  although  the  latter  tends  to  enhance  it  Figure  3c  shows  total  PV 
values  along  the  12.9*  C  isotherm  (the  200  dbar  temperature  at  y  -  60km),  which  passes  through  the 
minimum.  At  59*  W,  this  330  m  axis  minimum  is  absent  in  the  total  PV,  though  in  fact  the  sum  of  the 
stretching  and  non-QG  terms  alone  displays  it;  but  there  is  a  deeper  minimum  centered  at  700  m  about 
30  km  south  of  die  axis,  a  feature  hinted  at  in  Fig.  3a  for  63*W  by  the  slight  pinching  of  PV  contours 
at  that  same  depth,  again  30  km  offshore  of  the  axis.  Below  temperatures  of  7°-8*C  at  both  sections, 
PV  remains  approximately  uniform  along  isotherms. 

First  consider  how  well  the  model  reproduces  observed  PV  in  the  GS.  In  Figure  7,  the  PV  sec¬ 
tion  of  HF93  at  68*W  has  been  replotted  and  smoothed  using  three  contours  (10,  20,  and 
23  x  10~u  m~1a~1)  to  highlight  its  key  features;  it  is  shown  over  a  (yx)  range  comparable  to  that  in 
Figures  3  and  6.  (The  current  axis  y  -0  corresponds  to  the  location  of  station  7  in  HF93’s  figures.) 
Although  the  plot  is  noisy,  there  is  a  dear  break  in  strong  values  near  the  axis,  though  it  is  some  100  m 
deeper  than  the  corresponding  minimum  in  the  synthesized  section  at  63°  W.  HF93  identified  this  sin¬ 
gle  minimum,  and  associated  it  with  densities  of  "26-4-27.1  kg  m~3  (0  -  11"-18°C).  Figure  7  of 
this  note,  rather  less  cluttered  than  Fig.  9a  of  HF93,  suggests  that  the  minimum  occurs  in  a  deeper 
range  of  o0  (perhaps  26.9-27.4  kg  m~3, 6  «  8° -13.3°  C);  moreover,  at  40-30  km  offshore  of  the  axis 
between  600  and  700  m,  there  indeed  appears  to  be  another  relative  minimum  like  the  one  seen  in  the 
model  PV  at  39°  W. 

Relative  to  die  density  range  and  current  axis,  the  location  of  the  shallower  minimum  also  agrees 
dosely  with  that  pictured  in  Bower  et  al.  (1983).  Their  PV  section  at  68°  W  (their  Figure  Id),  calcu¬ 
lated  only  from  density  data,  shows  the  axis  minimum  and  offshore  relative  maximum  at  densities 
ofe  ■  27.0  -  27.3  kg  m~3.  They  do  not  find  die  weaker  and  deeper  secondary  minimum,  but  this  is  not 
surprising,  as  this  feature  depends  on  the  inclusion  of  terms  other  than  die  stretching  (compare  Figure 
6a  and  6b).  For  the  modd  at  63*W,  the  PV  contrast  between  the  minimum  and  the  higher  values  to 
north  and  south  (along  an  isotherm)  is  approximately  equal,  about  4-3  x  I0~nm~ls~l  (Fig.  3c).  The 
observations  at  6S*W  show  the  contrast  going  onshore  (northward)  to  be  much  greater  than  that  going 
offshore,  about  23  x  10~llm~ls~l  compared  to  3-10  x  1 0*lIm”,s“l  (HF93).  The  discrepancy  of  course 
is  due  to  die  feet  that  die  model  does  not  reach  the  high  values  of  PV  northward  of  the  axis,  as  noted 
above.  It  is  further  interesting  that  according  to  observations  die  shallow  minimum  does  not  seem  to 
persist  downstream:  it  is  barely  evident  in  Bower  et  al.'z  64"  W  section,  and  absent  entirely  in  their 
sections  farther  downstream;  nor  was  it  found  in  the  35°  W  CTD  section  (see  HF93).  This  is  consistent 
with  its  absence  in  Hendry’s  model  PV  structure  for  39" W  (Fig.  6a),  though  the  weaker  offshore 
minimum  at  y  s  30-60  km  persists.  As  noted  above,  the  latter  is  due  to  the  combined  effect  of  the  last 
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two  *■»—  on  tbe  right  hand  ode  of  (7),  both  of  which  are  due  to  the  flow  and  are  negative  south  of  the 
axis. 

We  return  to  the  dynamical  importance  of  this  specific  PV  structure  momentarily,  but  consider 
first  how  errors  in  the  fits  could  modify  die  results.  The  similarity  of  the  parameters  in  the  fits  for  A, 
D,  and  C  at  the  two  locations  reflect  the  persistent  downstream  coherence  of  the  GS’s  thermocline 
structure.  Indeed,  one  could  do  a  reasonable  job  of  selecting  these  fits  simply  by  careful  examination 
by  eye  of  the  observed  temperature  structure.  Quantitatively,  as  well,  the  model  is  insensitive  to 
changes  in  these  fits:  substituting  the  59°  W  fits  for  A,  D,  C  st  65®W,  the  rms  error  over  the  36  sta¬ 
tions,  between  200  and  1200  dbar,  increases  only  to  0-59°  C  (from  0.56°  Q.  Qualitatively,  a  substitution 
for  A  and  D,  or  for  C,  alone  does  not  affect  the  PV  structure  substantially. 

By  contrast,  changes  to  the  fit  for  file  local  inverse  scale  depth  B  more  strongly  affect  the  model. 
Quantitatively,  holding  the  fits  for  A,  D,  C  but  substituting  59°  W  parameters  for  B  increases  die  rms 
error  in  temperature  to  0.62°  C.  More  significantly,  this  substitution  completely  eliminates  the  relative 
minimum  at  330  m  seen  in  the  model  PV.  To  examine  this  point  further,  the  constant  parameters  in  the 
fit  for  B  were  varied  smoothly  so  that  the  steeper  65°  W  curve  in  Figure  4  gradually  approached  the 
shape  at  39°  W,  while  holding  the  fits  for  A,  D,  and  C  at  their  65°  W  values.  The  relative  minimum  in 
PV  persists,  retreating  upward  and  northward  as  the  59®  W  curve  for  B  is  approached  However,  the 
values  of  (a,  b,  c ,  d)  in  (2)  must  be  changed  over  90%  of  the  way  from  the  65°  W  to  the  59°  W  values 
before  the  minimum  disappears  altogether.  Indeed,  file  fitted  curve  for  B  as  given  by  the  parameter 
values  in  Table  1  is  significantly  better  than  a  fit  using  altered  values:  when  altered  values  (90%  of  the 
way  towards  the  59®  W  values)  are  used,  the  standard  deviation  of  the  difference  between  predicted  and 
fitted  values  of  B(Q  over  31  points  is  50%  greater  than  when  the  values  in  Table  1  are  used  Further¬ 
more,  over  most  of  the  range  of  C,  the  two  curves  differ  by  an  amount  greater  than  the  standard  devia¬ 
tion  of  die  original  fit 

We  now  return  to  file  dynamical  significance  of  features  in  the  PV  field  It  was  noted  that  the 
current’s  stability  properties  depend  on  its  PV  structure.  In  particular,  away  from  the  influence  of  top 

and  bottom  boundaries,  the  PV  gradient  45.  must  change  sign  somewhere  within  the  current  for  insta- 

dy 

bilities  to  be  possible.  Now  in  a  flow  where  density  is  conserved  following  fluid  parcels,  the  appropri¬ 
ate  gradient  II,  to  consider  is  that  calculated  along  isopycnals  (in  quasi geostrophic  theory  this  is  the 
same  to  lowest  order  as  n,  along  constant  z,  but  in  the  GS  isopycnals  depart  from  horizontal).  Follow¬ 
ing  file  PV  structure  along  the  isotherms  superposed  in  Figure  5a,  we  see  that  at  temperatures  below 
about  9®C  and  above  14®C  IT  is  nearly  uniform  along  isotherms.  However,  mid-thermocline  isotherms. 
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such  as  that  shown  in  Fig.  5c,  pass  through  some  part  of  die  minimum,  and  hence  along  these,  the  gra¬ 
dient  of  n  does  change  sign 

Theoretical  work  by  Pratt  et  al.  has  investigated  the  dynamics  of  a  current  with  such  a  PV  struc¬ 
ture,  modeling  the  flow  as  a  mid-latitude  jet  with  piecewise  uniform  PV.  [Note  that  because  the  axis 
minimum  j*  present  in  the  stretching  component  alone,  it  is  a  fundamental  feature  of  the  GS  and  it  is 
thus  appropriately  included  at  lowest  order  in  such  a  model.]  They  found  that  when  the  PV  contrast  at 
the  southern  edge  of  die  minimum  is  0(20-50%)  the  magnitude  of  the  contrast  north  of  the  minimum, 
instabilities  develop  an  the  southern  side  of  the  current,  which  resemble  the  "warm  outbreaks"  described 
by  Comillon  et  al.  (1986),  in  which  warm,  shallow  Gulf  Stream  water  (as  opposed  to  slope  water)  is 
ejected  into  the  Sargasso  Sea.  These  outbreaks  ate  observed  only  south  of  the  stream,  and  occur  only 
west  of  65®W.  At  68°  W,  the  PV  contrast  falls  within  the  range  necessary  for  instability  (HF93), 
though  in  the  synthesized  section  at  65°  W  the  northern  contrast  is  underestimated  for  reasons  cited 
above.  At  59°W,  the  shallow  minimum  is  absent,  and  indeed  warm  outbreaks  do  not  seem  to  occur 
this  far  east  Indeed,  Pratt  et  al.  ’s  (1991)  model  is  inappropriate  for  investigating  the  dynamics  of  the 
deeper  minimum  found  in  the  synthesized  PV  at  59°W,  because  that  minimum  arises  from  terms  not 
accounted  for  in  the  model. 

4.  Conclusions  and  Implications 

Following  a  method  introduced  by  Hendry  (1988),  36  XBT  temperature  profiles  near  65°W  have 
been  used  to  model  the  GS  thermodine  there  between  200  and  1200  dbar.  The  model  allows  the  tem¬ 
perature  range  and  scale  depth  of  the  thermodine  to  vary  as  a  function  of  thermodine  depth,  which  in 
turn  is  a  function  of  cross-stream  position  in  the  current  This  65°W  model  section  has  been  compared 
with  Hendry’s  analogous  result  at  59®  W:  the  quantitative  similarity  of  most  of  the  constant  parameters 
in  the  fits  for  the  two  models  reflects  die  downstream  persistence  of  die  GS’s  bulk  canonical  structure, 
inducting  mininwnn  and  nwHmum  temperatures,  minimum  and  maximum  thermodine  depths,  predicted 
baroclinic  transport  relative  to  1200  dbar,  and  the  current’s  deformation  radius.  There  is  a  marked 
difference,  however,  in  the  cross-stream  variation  of  the  scale  depth  for  temperature  at  the  two  sections, 
which  leads  to  a  significant  difference  in  the  potential  vortirity  structure  computed  from  the  model  at 
die  two  locations.  Ia  particular,  at  65°  W  there  is  a  relative  minimum  along  mid-thermocline  isopycnals 
near  the  GS  axis  at  350  dbar,  that  does  not  appear  in  the  downstream  section.  This  minimum  is  evi¬ 
dently  a  robust  feature  of  the  GS’s  PV  near  68°W,  inherent  in  the  density  structure  of  the  current,  and 
enhanced  by  the  inclusion  of  a  ncn-quasigeostrophic  term  in  the  expression  for  PV.  Overall,  the 
modeled  PV  structure  for  die  GS  at  65°  and  59°W,  and  the  downstream  modification  they  imply,  are 


conmtwtt  with  observed  PV  sections  between  68*  and  55*W.  At  59°W  there  is  a  deeper  relative 
minimum  30  km  south  of  the  axis,  not  seen  in  the  observations  of  Bower  et  al.  (based  on  density 
alone),  because  it  derives  from  die  relative  vortidty  and  noc-quasigesotrophic  contributions  to  the  PV. 

The  difference  in  the  types  and  locations  of  relative  minima  and  maxima  at  63*  and  59°  W  sug¬ 
gests  that  the  nature  of  instabilites  arising  in  the  GS  can  vary  as  a  function  of  longitude.  One  example 
is  the  warm  outbreaks  (Comillon  et  aL ,  1986)  observed  primarily  west  of  63°  W,  which  according  to 
theoretical  work  by  Pratt  et  al.  (1991)  could  arise  in  a  jet  with  the  PV  structure  observed  at  68*  and 
65°  W.  Alternatively,  the  cross-stream  variation  of  the  temperature  scale  depth  may  change  in  time  as  a 
result  of  meanders  or  ring  interactions  distorting  the  thermocline,  steepening  it  or  spreading  it  out 
Manning  and  Watts  (1989)  have  identified  two  important  modes  of  thermocline  variability  other  than 
translation,  that  could  have  such  an  effect  Their  so-called  "transport  mode”  is  associated  with  a  net 
increase  in  the  drop  of  the  thermocline  across  the  current,  and  their  "vorticity  mode”  with  a  change  in 
horizontal  shear.  The  first  would  certainly  change  the  fit  for  C(y),  and  the  second  might  affect  any  of 
the  fits.  Such  changes  are  likely  to  affect  the  axis  minimum  in  PV,  suggesting  that  temporal  variability 
in  die  current’s  stability  properties  could  result  as  well 
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Figure  1.  ft)  Harvard  Gulfcast  for  March  23,  1988,  with  EN175  cruise  track  overlaid 

b)  Temperature  profiles  for  the  36  XBT  drops  used  in  die  analysis,  c)  Temperature  contoured 
along  the  line  of  36  XBT  drops  shown  in  a).  Northeast  is  to  the  left,  southwest  to  the  right; 
cross-stream  distance  is  equal  to  distance  along  the  ship’s  track,  which  was  virtually  straight 

Figure  2.  a)  Ht  of  parameter  A  (eqa  (1»  to  C,  as  given  by  eqn  (2).  Pluses  mark  the  values  individu¬ 
ally  derived  for  each  of  the  31  stations  used  in  the  fit  Open  circles  are  the  values  derived  for  the 
5  stations  omitted  in  the  fit  b)  As  in  a),  but  a  fit  of  parameter  B  :o  C,  as  given  by  eqn  (4).  c) 
Fit  of  parameter  C  to  y  in  eqn.  (3).  d)  As  in  a),  but  a  fit  of  parameter  D  to  A,  as  in  eqn  (3). 

Figure  3.  a)  Temperature  contoured  for  the  southernmost  XBT  section  in  Fig.  la,  where  y  =  0  is  the 
station  location  where  T  »  13 °C  at  200  dbars.  (Positive  y  values  are  onshore  and  are  shown  to 
the  left.)  b)  The  model  temperature  contours  over  a  similar  range  of  y. 

Figure  4.  Comparison  of  the  fits  at  63®  and  59®W  for  B,  shown  explicitly  as  a  function  of  the  cross¬ 
stream  position  y.  The  inverse  of  B,  which  is  the  local  scale  depth  for  temperature,  is  plotted 

Figure  3.  a)  Total  potential  vorridty  field  for  the  13  parameter  model  derived  in  the  text  for  65°  W. 
Units  are  10"llm“ls_l,  and  the  contour  interval  is  2  x  10 ~ltm~ls~l.  Shown  as  dashed  contours  are 

the  6°,  10®C,  13®  and  16®C  isotherms,  b)  Stretching  PV  component, 

Po  dST 

Figure  6.  Same  as  in  3,  but  for  Hendry’s  parametric  model  at  39®W. 

Figure  7.  Observed  PV  at  68°  W,  adapted  from  Hall  and  Fofcnoff  (1993),  and  displaying  only  the 
10, 20,  and  23  x  contours.  Solid  areas  are  greater  than  25  x  hashed 

areas  are  20  -  23  x  l(T,lm"1*"1;  dotted  contour  is  the  10  x  lff'nm”1s_1  contour.  Dashed  con¬ 
tours  are  isopycnals  from  o®  -  26.9  kg  m~3  to  oe  -  27.4  kg  m~3  with  a  contour  interval  of 
0.1  kg  m~3. 


65** 


so®  iv 


T 

10.84  *C 

8.92  *C 

5 

0.00958  dbers-1 

0.00912  dbars~l 

a 

7.326  *C 

1230  °C 

a 

0.938 

0.942 

P 

4.829  *C 

4.830  *C 

a 

0.00170  dbars  ~l 

0.00265  dbars  “* 

b 

0.00315  dbars  "l 

0.00251  dbars  “* 

c 

658  dbars 

557  dbars 

d 

000216  dbars  -* 

0.00098  dbars 

326  dbars 

322  dbars 

By 

-0.0260  ton  _1 

-0.0248  ton 

Cy 

-9.74  ton 

-10.60  ton 

*>y 

437.4  dbars 

445  dbars 

Table  L  Values  of  the  13  parameters  in  (2)  -  (5).  The  column  labelled  65**  is  for  the  XBT  data  dis¬ 
cussed  in  this  note;  59*  W  refers  to  Hendry’s  (1988)  parameters. 


Table  2.  Comparison  of  transports  in  depth  ranges,  rr  :  maximum  velocities,  for  the  three  XBT  cross¬ 
ings,  and  for  the  model  sections  at  65°  and  59°W. 
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